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Using tlie liigli temperature approximation we study, within the linear response theory, the wake 
in the quark-gluon plasma by a fast parton owing to dynamical screening in the space like region. 
When the parton moves with a speed less than the average speed of the plasmon, we find that the 
wake structure corresponds to a screening charge cloud traveling with the parton with one sign fiip 
in the induced charge density resulting in a Lennard- Jones type potential in the outward flow with a 
short range repulsive and a long range attractive part. On the other hand if the parton moves with 
I a speed higher than that of plasmon, the wake structure in the induced charge density is found to 

, have alternate sign flips and the wake potential in the outward flow oscillates analogous to Cerenkov 

' like wave generation with a Mach cone structure trailing the moving parton. The potential normal 

CN| ' to the motion of the parton indicates a transverse flow in the system. We also calculate the potential 

due to a color dipole and discuss consequences of possible new bound states and J/tp suppression 
O ' in the quark-gluon plasma. 

o 

m 
(N 

(N 
> 



cn 
o 
o 



PACS numbers: 12.38.Mh,24.85.-|-p 



I. INTRODUCTION 



A plasma is a statistical system of charged particles which move randomly, interact with themselves and respond to 
external disturbances. Therefore, it is capable of sustaining rich classes of physical phenomena. Screening of charges, 
damping of plasma modes, and plasma oscillations are important collective phenomena in plasma physics 0, Q- A 
proper description of such phenomena may be obtained if we know how a plasma will respond macroscopically to a 
given external disturbance. The microscopic features of the particle interactions in the plasma are not completely lost 
in such a macroscopic description. But they can be implemented in a response function through the ways in which 
Clni the mutually interacting particles adjust themselves to the external disturbance and the response function plays a 
^ • crucial role in determining the properties of the plasma P|. 

Soon after the discovery of Quantum Chromodynamics (QCD), it has been found that at high temperature T the 
color charge is screened [Sj and the corresponding phase of matter was named Quark-Gluon Plasma (QGP). It is a 
special kind of plasma in which the electric charges are replaced by the color charges of quarks and gluons, mediating 
the strong interaction among them. Such a state of matter is expected to exist at extreme temperatures, above 150 
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■ MeV, or densities, above about 10 times normal nuclear matter density. These conditions could be achieved in the 
early universe for the first few microseconds after the Big Bang or in the interior of neutron stars. In accelerator 
experiments high-energy nucleus-nucleus collisions are used to search for the QGP. These collisions create a hot and 
dense fireball, which might consist of a QGP in an early stage (less than about lOfm/c) 4j. Since the masses of the 
lightest quarks and of the actually massless gluons are much less than the temperature of the system, the QGP is 
an ultrarelativistic plasma. To achieve a theoretical understanding of the QGP, methods from quantum field theory 
(QCD) at finite temperature are adopted Perturbative QCD should work at high temperatures far above the 

phase transition where the interaction between the quarks and gluons becomes weak due to asymptotic freedom. An 
important quantity which can be derived in this way is the polarization tensor describing the behavior of interacting 
gluons in the QGP. The dielectric function is related to the polarization tensor and important properties of the QGP, 
such as the dispersion relation of collective plasma modes and their damping or the Debye screening of color charges 
in the QGP 0. 

Recent numerical lattice calculations have found that charmonium states remain bound at least up to T ^ 2Tc Q 
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and the behavior of temporal correlators in pseudoscalar and vector channels deviates significantly from the free 
behavior at T ~ STc 0. These analysis suggested that pseudoscalars and vectors resonances may exist above Tq. 
Also robust results from Au+Au at the BNL Relativistic Heavy Ion Collider (RHIC) have shown collective effects 
known as radial To] and elliptic 0, 0| flows, and a suppression of high-pj^ hadron spectra Ts'l, which could 
possibly indicate the quenching of light quark and gluon jets (14j . The hydrodynamical description 15] of the observed 
collective flow indicates that the matter produced at RHIC behaves like a near-perfect fluid. On the other hand, the 
amount of jet quenching might depend on the state of matter of the flreball, i.e., QGP or a hot hadron gas. There 
are extensive theoretical efforts in understanding the effect of the medium on jet quenching 16, 17, 18, 19, 20, 21, 2^ 
using coUisional as well as radiative energy loss since the high energy partons traveling through a medium will lose 
energy owing to the interactions in the medium. 

Recently Ruppert and Miiller lE^I have studied the wake of the in-medium jet physics assuming a fast moving color 
charge particle with high momentum in a QGP within the framework of the linear response theory, considering two 
different scenarios: (i) a weakly coupled QGP as described by Hard Thermal Loop (HTL) perturbation theory and (ii) a 
QGP with the properties of a quantum liquid. In both cases the wake of the medium, i.e., QGP, has been observed 
through the induced charge and current densities due to dynamical screening but only in the latter case (ii) the wake 
exhibits a oscillation analogous to Cerenkov like radiation with a Mach cone structure. We point out here that their 
analysis within the HTL perturbation theory is not complete. In this paper we revisit their analysis within the HTL 
perturbation theory extending our investigations of dynamical screening |24j. 

In Rcf . 24] we restricted ourselves to the outward flow of the wake potential and to small particle velocities {v < 0.8). 
Also in the two-body potential for describing parton-parton interaction and bound states we did not include the 
magnetic interaction. Here we will show that both effects should not be neglected in the case of high-energy partons. 
Furthermore, in our earlier study we concentrated only on the screening aspect of the wake potential, whereas here 
we will discuss also the relationship to the energy loss and dynamic aspects (oscillations, flow). 

In particular, we flnd oscillations in the wake structure behind a fast moving color charge akin to Cerenkov like 
radiation and Mach cone. We further show that the wake potential for incoming flow is much like a Coulomb potential, 
while the potential in the outward flow flips its sign due to dynamical screening. In addition a potential perpendicular 
to the motion of the parton is also calculated which indicates a transverse flow in the QGP. We further compute a 
color dipole potential considering two moving charges and discuss the consequences of possible new bound states, viz., 
colored bound states, and J/tp suppression in the QGP. 



II. RESPONSE OF THE QUARK-GLUON PLASMA 

An appropriate description P, 0] of various plasma properties can be obtained if one knows how a plasma responds 
to an external disturbance. In order to establish a response relationship in a plasma, one usually considers the plasma 
response to the external electric field which induces a current density. If the system is stable against such a disturbance 
whose strength is weak, then the induced current may appropriately be expressed by that part of the response which 
is linear in the externally disturbing field. The linear response of a plasma to an external electromagnetic field has 
extensively been studied P, in plasma physics in which the external current is related to the total electric field by 

Kod^X). (1) 

where e{uj, k) is the dielectric tensor, describing the linear (chromo)electromagnetic properties of the medium and 1 
is the identity operator. Since we are interested in a system of a relativistic color charge moving through a QCD 
plasma, we apply the linear response theory for the purpose in the same spirit as in ^23] by simply assigning a color 
index, a = 1 • • • 8, to the relevant quantities here and also in the subsequent analysis to take into account the quantum 
and non-Abelian effects. Non-Abclian effects (beyond the color factors, e.g., in the Debye mass) will be important at 
realistic temperatures. Unfortunately, they cannot be treated by the method used by Ruppert and Miiller |23j and 
are therefore beyond the scope of this work as we closely follow their work in our analysis. 



e(w,k) - 1 



A. Linear Response in a (Color) Plasma 

In this subsection we briefiy outline the theoretical description P, based on the linear response theory for the 
QGP as a finite, continuous, homogeneous and isotropic dielectric medium, which can be expressed by a dielectric 
tensor, e(tj,k), depending on the direction only through the momentum vector, k. One can also construct another 
set of tensors of rank two from the momentum vector, k which are the longitudinal projection tensor, and the 
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transverse projection tensor, , respectively, given as 



(2) 



with the general properties ofthe projection operators: {V^f = V^, {V'^f = V'^ , V^-V"^ = and {V^f + {V'^f = 1. 
Now the dielectric tensor can be written as a linear combination of these two mutually independent components as 



e(c^,k) = V'-eLiio^k) + r^eTiio,k) , 
= Vt,eL{to,k) + Vj^eriio^k) , 

where the longitudinal, and the transverse, ey, dielectric functions are given by 



(3) 



eT{uj,k) 



fc2 



Treiuj, k) — 6^(0;, k) 



(4) 



Now, we shall use the dielectric tensor in (O for a set of macroscopic equations, viz., Maxwell and continuity equations, 
in momentum space to obtain a relation il| between the total chromoelectric field and the external current as 



lOJ 



(5) 



If there is no external disturbance to the plasma, J^^j = 0, (0 reduces to 

k^' 



ii;r„t(c.,k) = o. 



(6) 



which has nontrivial solutions representing dispersion relations in a linear medium, if the determinant vanishes, i.e 



det 



= , 



(7) 



leading to the following equations defining the longitudinal and transverse modes 

elIuj, fc) = , 

^2 

€T{uJ,k) = — . 



(8) 



So, the dielectric tensor contains essentially all the information of the chromoclectromagnctic properties ofthe plasma. 
However, one can also study a density response 1] to an external test-charge field through a density fluctuation which 
could induce a color charge density as 



£(w,k) 



1 P?xt(^,fc) , 



(9) 



where p^^t is the external charge density and e{uj,'k) is an another dielectric response function that can be obtained 
as 



k • e{uj, k) • k 



(10) 



This is analogous to the longitudinal dielectric function in I^J with the difference that it is independent of the direction 
of k because of isotropy. For a homogeneous and isotropic system, one can write 



eLiuj,k) = £(a;,k) , 



(11) 
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which can lead to a longitudinal dispersion relation for density fluctuations indicating that the space charge field could 
be spontaneously excited without any external disturbance. 

Combining Hll|l with ((HJ the induced color charge density by the external charge distribution becomes 



- 1 PcKti^^k) , 



and the total color charge density is given as 



eL(w,fc) 



(12) 



(13) 



in which the dielectric function measures the screening of an external charge due to the induced space charge in the 
plasma. 

The screening potential in momentum space can be obtained from Poisson equation in Coulomb gauge P, |3| a-s 



$"(w,fc) = An 



Pcxt{^,k) 



Combining ^ and jSJl the induced current density can be related to the external current density as 



^^d('^,k) = 



-i]r 



1 — er T 



tT — k"^ jup' 



^ext(^,k) 



(14) 



(15) 



B. A Fast Color Charge in a Plasma 



Since we are interested in studying the wake behavior of the plasma reacting to a charged particle, moving with 
a constant velocity v, we specify the external current and charge density, respectively, as 



J'l^^ = gVJ(x - vt) ^ 2^QVJ(co - k • v) , 



Poxt 



FT 



Q°5(x - v<) 27rg''(5(w - k • v) 



(16) 



where FT stands for Fourier transformation and v ~ |v|. The finite motion of the charged particle deforms the 
screening charge cloud and also suffers a retarding force causing energy-loss. The soft contribution to the coUisional 
differential energy-loss by the induced chromoelectric field is defined as 0, l25l | 



The induced electric field can be obtained from the inverse of (@J which is given as 



(17) 
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^eL J \eT — k'^l'-P l — k'^/u}'^ 

Combining and H18() the soft contribution to the differential coUisional energy-loss in QGP can be obtained |2 
as 

dE Ca 



(18) 



dx 



Ime^i + {v^k^ - uj^) Im {lj^^^ - k^' 



(19) 



where w = k • v originates from (|15|l . C is the quadratic Casimir invariant, and as — 5^/47r is the strong coupling 
constant. 



C. Response Functions in the High Temperature Limit 

The dielectric functions, both longitudinal and transverse, in are related @, H [H HI EHH to the self-energies 
of the gauge boson, i.e., gluon, in the medium as 

eLiu;,k) = 1 , 

nT('^, k) 

eT{u;,k) = 1- ' , (20) 
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where K = (w, k) with k = |k|. and IIt are, respectively, the longitudinal and the transverse self energies of the 
gluon, which are to be calculated in the high temperature limit. 

The transversality of the polarization tensor, K^Yi^n — 0, implies that only two components are independent, and 
the general solution can be written as 



n^^(w, k) = ^^^HtC^^, k) + B^^IVl{lo, k) 
where the transverse and longitudinal projection tensors [2^ are 



(21) 



1 



(22) 



in which is the four velocity of the fluid with u^u^ = 1. The tensor f)^^ and the vector K which are orthogonal to 



are defined as 



(23) 



with the Minkowski metric tensor rj^i,. Using the properties of the projection operators, one can obtain the scalar 
functions in H21|) as 



IlL{uJ,k) = -— noo(t^,fc). 



nT(a;,fc) = i(^5,^.-Mi^n,,(if) = i[n{;(c.,fc)-n4^,fc) 



(24) 



For computing the self energies one needs the in-medium gluon propagator which can be obtained from the Dyson- 
Schwinger equation 



(25) 



where AjJ^(i^) is the gluon propagator in vacuum. Combining H21|l . I|22|l and (|25|l . the full in- medium gluon propagator 
in covariant gauge ^ is obtained as 



■Ani/ 



if 2 -Ut k^- Hl 



(26) 



where ^ is the gauge parameter. 

The one loop gluon self energy in the high temperature limit has been obtained in Refs. [2^l30j which was found to 
be equivalent ^6, 31^ to the HTL approximation 32] that relies on the restriction to hard loop momenta and is given 
as 



ULiuj,k) - 

UT{uJ,k) = 

where the Debye screening mass follows from 



fc2 



1 



LU LU + k 

— In 

2k LU — k 



"2~fc2 



bJ ^ bJ + k 



k^ J 2k 



In 



mjj = nL(w = 0, fc) = noo(w = 0, fc) g^T^ 1 
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(27) 



(28) 



in which Nf is the number of light quark flavors in the QGP. Though the gluon propagator is a gauge dependent 
quantity, the gluon self energies in H27|) are gauge independent only for the leading term of the high temperature 



^ This particular gauge is chosen for convenience as we will see later in subsec. Ill Fl 
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FIG. 1: Dispersion relations for the longitudinal (L) and the transverse (T) part of giuons in the QGP (quadratic scale). The 
dashed line indicates the light cone, uj = k . 



expansion 
as 



The dielectric functions are therefore also gauge invariant and can be written combining (|2()|l and H27|l 

-k 



eL{uj,k) = 1 



k^ 



1 



2k 



In 



- k 



2 r, ,2 



eri.^k) = 1-^ 



k 



inQ{k-LU^) 



(29) 



The dispersion of longitudinal and transverse giuons, determined by the respective real parts of the dielectric 
functions in © (also can be obtained from the poles of the gluon propagator in H26|l ). are shown in Fig. ^ The 
following important features of the plasma are transparent from Fig. ^ (i) There are propagating modes above a 
common plasma frequency cu > ujp\ = yml^/S, which are always above the light cone {cu > k), show collective 
behavior, and are called normal modes. The longitudinal mode (plasmon) propagates with an average phase velocity 
■y/3/5c whereas the transverse mode with -^/G/Sc for moderate k [2^. (ii) In the static limit a; = 0, the self energies 
11^(0, k) — mj^ and IIt{0, k) = 0, respectively. This means that the longitudinal branch indicates the static screening 
in the electric sector with an inverse screening length, rrio, given in H28|l . whereas the solution of the transverse branch 
at A: = 0, implies the absence of magnetostatic screening, (iii) For < w < Wpi, k becomes imaginary |2^ and there 
are no stable normal modes as in (i), however, there is a collective behavior which corresponds to dynamical screening 
both in electric and magnetic sector as long as w > 0. (iv) For u; < k, the imaginary parts of the dielectric functions 
corresponds to Landau damping causing energy dissipation in the plasma through elastic collisions. 

When a color charge moves with a velocity, v, the dispersion relation disappears from above the light cone and 
owing to the interaction of the moving test charge with the medium, it is possible to define another mode as 



UJ 



kv 



(30) 



which has the support from due to the (5-function restriction. The collective modes in H30(l are in the space like 
region in which Re(ci;) < k for < w < 1 in contrast to the normal collective modes, where Re(w) > k always holds as 
discussed above in case (i). The dielectric response functions in H29|) provide a direct measure of the extent of which 
an external test charge is screened by the induced space charge in the plasma. The screening in this case is dynamical 
in the sense that it depends on the frequency as well as the wave vector similar to the case (iii) with the difference 
that the wave vector is damped |2^ there. 

Equation (|30|l is also known as Cerenkov condition 1] for wave emission by a moving parton in a plasma. Since 
the average phase velocity of the transverse modes yJ&JZc is greater than the speed of light for the isotropic plasma, 
the Cerenkov condition is not satisfied. According to (|3U|I there could be two important effects due to the particle 
interaction with the plasmon wave, which is propagating with an average phase velocity -y/S/Sc = 0.77c: 

1) The modes which are moving with a speed less than the average speed of the plasmon mode (0.77c) can be 
excited but they accelerate the slower moving charge particles and decelerate those moving faster than the modes. In 
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the former case energy of the excited modes is absorbed whereas in latter case energy is transfered to those modes 
by the charges. Such modes are excited by a charged particle which moves less than 0.77c but are severely Landau 
damped owing to absorption and emission resulting in a wake in the induced charge density as well as in the potential. 

2) The modes which are moving with a speed greater than 0.77c can be excited and they may not be damped, 
which are different from the existent plasma waves as discussed in case 1). When a charged particle moves faster than 
0.77c, it can excite undamped modes which could generate Cerenkov like radiation and Mach shock waves leading to 
oscillations in the induced charge density and in the corresponding potential. We note that Mach waves are generated 
in the plasma in general if the particle moves faster than the speed of sound, Cg in the plasma and the collective modes 
shift to the space like region. The Mach opening angle is given by 

(l)M — arccos • (31) 

In the following we would like to investigate those effects of the dynamical screening due to the fast moving charged 
particle on the wake of a color plasma in the high temperature limit. 



D. The Wake in the Induced Charge Density 



Combining H12(l and (|16|l . the induced charge density in configuration space reads as ^4 



In cylindrical coordinates k — 
density in 1(5^ becomes 



(fik 



dhJ 



(27r)3 

(kcos(/), Ksin0, fc^), f 



i • r — bjt^ 



(5(cj — k • v) 



(32) 



(p, 0, z) and assuming v — (0, 0, w), the induced charge 



n dn Jq{kp) 



+ C30 



duj exp 



(- 



1 



eL(w,fc) 



(33) 



where Jq is the Bessel function and k — + lu^ /v^. 
Using the symmetry properties of 6^(0;, k), viz., 



ReeL(— a;,fc) = RceL(ijj,fc), 
ImeL(-a;,fc) -ImeL(a;,fc), 



the induced charge density in (|33|l reads as 



(27r)2w 



K dn Jo (up) 



duj 



cos 



("(^-0) 



[ Re eL 
V A 



- 1 



sin I a; I t 

■.V 



Im eL 



(34) 



(35) 



where A — [Re cl]'^ + [Im e^]^ 

Now we would like to point out that, if the external charge is at rest, v — > 0, relative to the background, the induced 
charge cloud remains spherically symmetric. In the limit v — > 0, (|35|l reduces to 



47r 



Jq{kp) 



An 



: exp 



(36) 



This can be understood as follows: the background particles are moving isotropically on the average around the test 
charge. Introduction of a test charge which is at rest merely introduces a local fiuctuation of the number density in 
its vicinity but does not spoil the spherical symmetry of the system. As a result the response of the medium modifies 
the induced charge density into a static Yukawa type still refiecting this symmetry. 

The situation, however, changes if the test charge is in motion relative to the heat bath. The motion of the particle 
fixes the direction in space and spherical symmetry of the problem reduces to axial symmetry. This implies the loss 
of spherical symmetry of the Debye screening cloud around the moving test charge resulting in a wake in the induced 
charge due to dynamical screening as given in (|35|l . 

In Fig.|21we display the equicharge lines in the induced charge density for v = 0.55c in cylindrical coordinates which 
are scaled with rrig = \J (The induced charge density is proportional to m^.) The induced charge density 
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FIG. 3: Left panel: The plot reproduces the equicharge lines in the induced charge density for v = 0.99c with a fixed color 
charge Q" in the high temperature plasma as obtained in Ref. [2^. Right panel: Same as left panel but plotted with more 
number of equicharge lines with enlarged axes to show the oscillations in the same density wake. 



is found to carry a screening charge cloud with a sign flip along the direction of the moving charge. As discussed 
in the preceding subsection III CI there wiU be an excitation of modes during such a wake but they will be severely 
Landau damped through induced emission and absorption of the plasmons by scattering off single particles. As a 
result the particle suffers an energy loss through elastic collisions which can be described by (|19|l and has been studied 
extensively in the literature 0|. When a charged particle moves with a velocity of less than 0.77c, < 2.5, the soft 
part of the colhsional energy loss, caused by the induced chromoelectric field, turns out to be important |2lll2a . l34Ll35l| . 

Next we consider the case when a charged particle moves faster than the average speed of the plasmon. For this 
purpose we choose v — 0.99c and the result is shown in the left panel of Fig. 13 which reproduces the corresponding 
results of Ref. |123| . The wake in the induced charge density apparently does not show any sign flip as in Fig. El and 
it shows merely a screenin g cl oud around a test charge similar to that given in H36|l . However, this is not the case as 
the particular analysis in |23l | was not complete and some important features were obviously missing in the wake of 
the medium. 

The complete analysis is shown in the right panel of Fig. |21 when the same data corresponding to the left panel are 
plotted with enlarged axes with a larger number of equicharge lines in the induced charge density. As expected the 
wake of the medium is clearly visible in the induced charge density which follows the charged particle traveling with 
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a velocity, v = 0.99c. Moreover, the screening cloud following the charge particle is found to flip its sign alternately 
indicating enhanced and depleted charge density. As discussed in the preceding subsection 111 CI there will be an 
excitation in the plasmon modes during such a wake which could emit a Cerenkov like radiation with a Mach cone 
structure trailing the moving parton. The radiation is like a spontaneous emission of plasmon waves and cannot be 
calculated through the expression given in H19|l and requires a special treatment. It is also important to note that for a 
parton velocity v — 0.99c, > 4, corresponding to the ultra-relativistic limit for which the hard contribution 

of the coUisional energy loss due to elastic scattering and the radiative energy loss are expected to be dominant. 

When a parton moves supersonically, it excites waves in a colored plasma. The propagation of sound waves could be 
estimated from the emission pattern of secondary particles traveling at an angle with respect to the jet axis through 
the relation pLIs^ 

A(/) = 7r±(/)M = TT ± arccos ^— ^ , (37) 

where = tt corresponds to the location of a maximum distribution of secondary hadrons from the away side jet 
m p + p collisions, where no medium effects are present. However, due to medium effects the azimuthal angle for 
secondaries is expected to be peaked differently from tt and could be related to the Mach opening angle, 0m- 

There are several programs in STAR and PHENIX experiments at RHIC where various approaches are being 
followed in order to find Mach cones. The basic approach is to use correlations between high p± secondary particles 
in pseudorapidity and azimuthal space. Some of the background correlation effects which might mimic Mach cones 
are being studied in details. Both the experiments are looking at 3-particle correlations which will eliminate the 
trivial backgrounds. Now already some interesting results have been reported by the PHENIX collaboration j^^l 
where the data suggest that the peak in the secondary particle correlation provides compelling evidence for a strong 
modification of the away-side jet. This could indicate two possible scenarios: (i) a Cerenkov jet in which the leading 
and away-side jet axes are aligned but fragmentation is confined to a very thin hollow cone or (ii) a directed jet in 
which the away-side jet axes are misaligned. A more quantitative investigation is necessary to distinguish between a 
'Cerenkov jet' that leads to a Mach cone, and a 'directed jet' which is misaligned. Following (|37|l the azimuthal 
angle can be estimated here as A0 « tt ± 1.01 using = 0.48 within the HTL approximation. Once a characteristic 
correlation structure in secondary hadrons leading to a Mach cone, 4>m ~ 1-01 is confirmed experimentally, it would 
provide an estimation of the angular structure of the energy loss and also the speed of sound in the QGP. 

In Ref.li^ it was argued that the wake of an energetic parton can create an observable flow in the QGP. However, 
it was also argued that this is only the case if the energy loss of the parton is very large (about 12 Gev/fm). Even 
combining the coUisional and radiative contributions, we do not expect such an high energy loss in relativistic heavy-ion 
collisions. 



E. The Wake in the Screening Potential 



Combining H14|) and H16|) the screening potential in configuration space due to the motion of a color charge can be 
written as ,24] 



j,ik-(r— vt) 



In cylindrical coordinates, the screening potential in (|38|l becomes 



(38) 



(p,z,i) 



2Q" 



dn kJq {up) 



dco 



1 



fc2A(w,fc) 



COS |a> (~ ^ ^) } ^ ^^^^ (~ ^ ^) } ' ^"^^^ 



where, k — ^ + uj^ /v^ and A = (Rcel)^ + (Ime^)^. 

The potential in (|38|l is also solved for two special cases, (a) along the direction of motion of the parton (r || if) and 
(b) perpendicular to the direction of motion of the parton (r _L w). The potential for the parallel case is obtained as 



$ff (r, V, t) 



2g" 



oo /•! 

dk I dx 
Jo 



cos [kx |r — vt|) — h sm [kx |r — vt|) — - — 



(40) 



whereas that for the perpendicular case is 
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FIG. 4: Left panel: Spatial distribution of the scaled wake potential with respect to mo from a jet with a fixed color charge 
Q" which is traveling with v = 0.55c. Right panel: This plot shows the corresponding equipotential lines. 
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FIG. 5: Same as Fig. H but v = 0.99c. 



Qa roc i-^TV 

$1 (f, v,t) ^ ^ I dk I d(j) I dx 



2tz /•! ifc r cos (p^/l — x'-^ — vxt 



27r^ Jo Jo J-i (t^ vkx, k) 



/ dk dx — [cos (kvxt) Re — sin [kvxt) Im e^^] , (41) 

Jo Jo ^ 



with X — cos 6*. We point out here that some efforts were made earlier in studying the wake potential [2J, |39j in a 
QGP. However, those were incomplete, neglecting Mach cones and Cerenkov radiation in the wake of the plasma due 
to the motion of the charged particle. 

In Fig.^lthe wake potential is plotted for v = 0.55c in cylindrical coordinates which is scaled with mo- (The wake 
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FIG. 6: Left panel: Scaled screening potential along the jz-axis, i.e., parallel to the direction of the moving color charge for 
different velocities. Right panel: Same as left panel but along the p-axis, i.e., perpendicular to the direction of the color charge. 



potential is proportional to mr).) The left panel shows a three dimensional plot of the wake potential whereas the 
right panel displays the corresponding equipotential lines. It should be noted that the potential depends only on v 
and not on t as it should be for an isotropic and homogeneous plasma. The spatial distribution of the potential shows 
the well known singular behavior at r = 0, i.e. z = and p = 0, which has been cut by hand through the scale 
restriction in the left panel. It also exhibits a well defined negative minimum in the {p — z) plane which has been 
formed due to the dynamical screening in the vicinity of the moving charged particle. This feature is also obvious 
in the right panel where the equipotential lines are found to flip their sign in the (p — z) plane corresponding to a 
moving induced charge cloud of opposite sign as shown in Fig. |21 forming the wake around a moving particle. 

When the velocity of the charged particle is higher than the average plasmon speed (0.77c), the induced charged 
cloud trailing the charged particle begins to oscillate which in turn indicates an oscillatory wake potential and has not 
been found in earlier studies [23.l39l |. In Fig.|31we display the spatial distribution of the wake potential for v — 0.99c 
which shows oscillations in the (p — z) plane. Such oscillations, as discussed in the preceding sub-sections, lead to 
an emission of Cerenkov like radiation and generate Mach shock waves. This is analogous to the wake phenomena in 
classical plasmas in the supersonic regime pOj . 

In Fig. El the wake potentials in the specific directions corresponding to Figs. 0| and |31 are displayed. The left panel 
shows a wake potential along the 2;-axis, i.e., parallel to the direction of the moving color charge. This could be 
directly obtained from H40|l or by setting p = in H39|) . As expected the singular behavior at z = is clearly reflected 
in this plot. In the outward flow, i.e., in the negative z-direction the wake potential compared to the static one (i.e., 
for V = 0) falls off very fast, reverses its sign, exhibits a negative minimum and asymptotically approaches zero from 
below. For v < 0.77c the potential in the outward flow is found to be of the Lennard- Jones type Md] which has a short 
range repulsive part as well as a long range attractive part. With increase of v the depth of the negative minimum 
increases and its position shifts towards the origin or to the particle. When 0.77c < v < c, the wake potential begins 
to oscillate. For v = 0.99c such oscillations in the wake potential along the direction of the motion are also clearly 
visible in the outward flow (z < 0) from the inset of the left panel in Fig. |S1 On the other hand the potential for 
incoming flow^, i.e., in the positive z-direction does not show any such structure and behaves more like a modified 
Coulomb potential. However, with the increase of v it attains a Coulombic form indicating that the forward part of 
the screening cloud is not so strongly affected by the motion of the particle. 

Figs. 01 13 and El reveal that at finite v the potential in the z-direction becomes anisotropic and loses forward- 
backward symmetry with respect to the motion of the charged particle. The origin of this can be traced back to the 
second term in (|39|l and/or l|40|) which is antisymmetric under inversion of [z — vt) —t—{z~ vt). This term appears 
due to the antisymmetric nature of the imaginary part of eL(k,uj) in H34|) . So, the effect of finite v is not limited only 
to the deformation of the screening cloud, but, more importantly, it shifts the position of the test charge off the centre 



^ Our earlier study |24| was restricted only to the outward flow and small particle velocities. 
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of the wake potential. As a consequence, the test charge experiences a retarding force from the plasma. This leads to 
a stopping power of the QGP against a moving test charge. The nature of the energy loss leading to such stopping 
power of the QGP has already been discussed in the preceding subsec. IIIDI We also note that a minimum in the 
screening potential for a moving color charge was found in earlier works by Chu and Matsui and also by us 
in the direction of propagation. However, in Ref. 39] Chu and Matsui did not report this minimum whereas our 
earlier study was incomplete concerning the forward-backward asymmetry of the wake structure and Ccrenkov 
radiation. 

On the other hand, the right panel of Fig. displays the wake potential normal to the direction of motion of 
the charged particle. This could be obtained directly from |(1T|) or by setting z = in which shows the usual 
singularity at p = and a symmetric behavior in p. It also falls off very fast compared to the static case and the 
nature of the potential is found to be of the Lennard-Jones type due to the deformed screening cloud forming the 
wake in the vicinity of the moving charged particle. With increase of v, the position of the minimum shifts towards 
the origin and the depth of it increases. Such kind of potential in the normal direction implies that a part of the 
screening cloud forming the wake due to the motion of the charged particle along the ^-direction is also moving away 
in the perpendicular direction, resulting in a transverse flow in the medium. This behavior was not observed in earlier 
studies .24,39,]. 

The negative minimum in the wake potential indicates an induced space charge density of opposite sign. Thus, a 
particle moving relative to a particle in the induced space charge density would constitute a dipole oriented along the 
direction of motion. In the next sub-section, we calculate the potential due to such dipole interaction in the QGP. 



F. Dipole Potential in the QGP 



We consider two color charges and separated bya distance r. The change in free energy of the system to 
bring the two widely separated color charges together [3, |43, 13 is 



where J^^^ {x) = 



Pcxt 1 " cxt 



duj 



2 J {2^Y J-oo (27r) y_eo (2^) 



(42) 



is the sum of the two external currents, (x) — Jf (x) + J2 (x), and (x) is the 

associated gauge field. The entropy generation is neglected in 142() . 

We assume that the average value of vanishes in equilibrium (^/j)cq 
vector potential, A^^ follows from linear response theory j2^ as 



0. The induced expectation value of the 



(43) 



where, Q = u.K, k = \/iW—K^ are the Lorentz invariant energy and 3-momentum, respectively. A^^ is the 
propagator of the gauge boson exchanged between the two currents which is given in 1)26(1 . Combining H43() with 142() 
we can write 



AT = 27r 



d^k 
(2^ 



duj' 



{2n) (27r) 



J^,t c.,-k A^,(c^',k)J,-,t(c^',k). 



(44) 



One can obtain the dipole interaction in H44|l by separating |2^ the current"^, in general, into a charge density that 
moves with the fluid velocity and a spacelike current flow either longitudinal or transverse to as defined in 
(113: 



(45) 



where 



K^'uo 



J\K). 



(46) 



^ Henceforth, wc drop the suffix 'ext' from the external current for convenience. 
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Using H26|l and (|45|l the dipole interaction can be obtained as 



JT(^,-k)-JT(t^,k) 

Louj' - fc2 _ni,(w',fc) - fc2 - nT(t^',fc) 

3{u;, -k)-J(co, k) + ^pju;, ~k)p{cj, k) 
u'^ -k^ -IlT{uj\k) 
k'^ + ujuj' — k)p(aj, k) 

k^ Lu'^ - fc2 _nL(a;',fc) ' 



(47) 



where u.J p. Though the gluon propagator in H2t)|l is gauge dependent, the dipole interaction becomes gauge 
independent. 

We begin with the known example of screening of two static charges and at position xi and X2 , respectively. 
The corresponding dipole current is given by 



(x) = Jf [x) + Jl^ {x) = 5^ [Q''S (x - xi) + Q^S (x - X2)] ^= 2nS^S (u;) Q'' 



-ikx2 



(48) 



Combining (|17|l and and also dropping the self coupling terms we get the well known Yukawa potential as 



27r2 



d-^k- 



(k.f) 



(49) 



where, r = xl — xj and r — (p,0, z). It is also noteworthy to mention that in Ref. f43| a static screening potential 
was also obtained in a QGP, where a polarization tensor beyond the high temperature limit was used. However, this 
approach in Ref. |42| has its limitation because a gauge dependent and incomplete approximation for the polarization 
tensor was used. 

Now for two comoving charges and the dipole current can be written as 



J^(t,x) = (l,v) [Q"(5 (x - xi - vi) + Q''<5 (x - X2 - vi)] 2^ (1, v) 5 (co - k • v 
and we obtain the dipole potential as 



— ik-X2 



(50) 



27r2 



AT'^^r-p, z) 



(Pk cos ^k • 

dn kJq {up) 



- Ut {lu, k) - Ul (w, k) 



dkz cos {kzz) < Re 



fe2 



1 - ^ 



- Ht {uj, k) - Ul (w, k) 



(51) 



The following features are transparent in H51|l . The first term corresponds to the transverse (magnetic) interaction 
whereas the second term is due to the longitudinal (electric) interaction^. Like the single body potential in (|39|) the 
two body potential is also symmetric in p. Moreover, in (|51|l both imaginary parts corresponding to the longitudinal 
and transverse response functions drop out because they are odd in kz , which make both terms symmetric under the 
inversion of z ^ — z, thus the total potential. This is, however, unlike the single body potential in H39|l where the 
presence of an imaginary part of the longitudinal response function was found to be responsible for the asymmetric 
behavior. In addition the two body potential should be symmetric in z as the two color charges are moving opposite 
to each other which is obvious from the dipole interaction in H47|) . Recalling Ampere's law this would also amount to 
an attraction between color charges traveling opposite to each other along the z axis. Both the electric and magnetic 
interactions play a crucial role, as we will see below. 

In Fig.|7|the spatial distribution of the scaled ^ longitudinal part of the two body potential in cylindrical coordinate 
is displayed for v = 0.55c. The left panel shows a three dimensional plot of it whereas the right panel displays the 



In our earlier study the two body potential was obtained by averaging the single body potential and found to depend only on the 
longitudinal (electric) interaction. 

^ The potential is scaled with mo as well as with the interaction strength, Q"^Q^. So, the details of the potential will depend on the 
temperature, the strong coupling constant, and the sign of the interaction strength. 
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FIG. 7: Left panel: Spatial distribution of the longitudinal part of the scaled dipole potential with respect to mu in which 
the two color charges are traveling antiparallel to each other with v = 0.55c. Right panel: This plot shows the corresponding 
equipotential lines. 




corresponding equipotential lines. The longitudinal part has a singularity at r = 0, i.e., z = and p = 0, clearly 
reflected in Fig. |7| As discussed above the potential plotted in both panels is found to be completely symmetric in 
the {p — z) plane along with a variable minimum. Such a minimum in the {p — z) plane is due to the fact that the 
electric interaction contributes differently in z as well as in p direction. The potential has a positive minimum in the 
longitudinal direction, i.e., along the z axis and then tends to zero for large \z\. This can also be seen from the left 
panel of Fig. 1101 which can be obtained from (|51() by setting p — Q. On the other hand, in the transverse direction, 
i.e., along the p direction, the electric contribution has a well defined negative minimum, it also attains a positive 
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FIG. 9: Total contribution for v = 0.55c 
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FIG. 10: Left panel: Scaled two body potential along the direction of motion when two color charges are traveling antiparallel 
to each other with v = 0.55c. Right panel: Same as left panel but normal to the direction of motion. 



maximum at large p and slowly approaches zero asymptotically. This feature can clearly be seen from the right panel 
of Fig. EH which is obtained by solving H51|l with z = Q. 

The contribution of the scaled transverse part is shown in Fig. |S1 which is also found to be symmetrical in the (p— z) 
plane along with a well defined negative minimum. However, both panels indicate that the magnetic interaction 
contributes differently in z as well in p directions, which can also be seen from Fig. 1101 In the limit u ^ 0, the 
transverse part vanishes. However, for any nonzero value of v the transverse part begins with a finite negative value 
for both cases. As \z\ increases, its magnitude gradually decreases and approaches zero for large \z\ (left panel of 
Fig. I10() whereas in the transverse direction its magnitude decreases faster, becomes positive at certain values of p, 
and then slowly approaches zero for large values of p (right panel of Fig. I10|l . 

The spatial distribution of the total potential for v = 0.55c is displayed in Fig. which appears due to the 
compensating effects between the electric and the magnetic interactions as discussed above. The resulting potential 
shows the usual singularity of the screening potential at r = (z = and p = 0) and a completely symmetric behavior 
along with a pronounced negative minimum in the {p — z) plane, which are clearly reflected in both panels of Fig. 
However, the detailed features in the specific direction can also be seen from both panels of Fig. ^| The potential 
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FIG. 11: Left panel: Spatial distribution of the longitudinal part of the scaled dipole potential with respect to mo in which 
two color charges are traveling antiparallel to each other with v = 0.99c. Right panel: This plot shows the corresponding 
equipotential lines. 




-10-5 J LO 



FIG. 12: Same as Fig. Illl but showing the transverse part. 

along the |z|-direction as shown in the left panel falls off like that of a static one, flips its sign, exhibits a negative 
minimum, and then oscillates around zero. As can be seen from the left panel, the competing effects between electric 
and magnetic contributions with opposite sign result in an oscillatory potential along the z direction. On the other 
hand, the potential in the transverse direction (right panel) falls off slowly compared to the static one, flips its sign, 
and exhibits a negative minimum at some values of p, again changes its sign to attain small positive maximum and 
then tends to zero at large p. It can be seen that the nature of the potential is completely dictated by the electric 
interaction rather than by the magnetic one. In both directions there are oscillations at large distances but the form 
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FIG. 13: Total contribution for v = 0.99c 




FIG. 14: Same as Fig.[T01but i; = 0.99c. 



of the potential mostly resembles the Lennard- Jones type with a pronounced repulsive as well attractive part. 

With the increase of v, as we will see below, the relative strength of the dipole potential grows strongly. For 
V = 0.99c Fig. ^2 displays the spatial distribution of the scaled longitudinal potential whereas Fig. 1121 shows that of 
transverse part in the {p — z) plane. These plots clearly show that the magnitude of the spatial distribution of the 
dipole potential due to both electric and magnetic interactions are enhanced to a large extent in the (p — z) plane. 
The resulting potential distribution is shown in Fig. ^] Apart from the usual features as discussed above for smaller 
w, it displays a new and important feature. A substantial repulsive interaction has grown in the transverse plane, i.e., 
in the direction of p but at z = 0, becoming responsible for a vertical split of the minimum in the (p — z) plane as 
can be seen from both the panels in Fig. E| 

Now, the detailed contributions of both the interactions can also be understood from both panels of Fig. ^] where 
the dipole potentials for two special cases have been displayed. As shown in the left panel of Fig. E|the depth of the 
negative minimum of the dipole potential along the \z\ direction increases and it position shifts towards the centre 
indicating a faster fall off as well as a larger anisotropy. This is due to the fact that the magnitude of the magnetic 
interaction becomes dominant. The form of the potential remains of the Lennard- Jones type. On the other hand, the 
potential in the transverse direction (z = 0) as shown in the right panel of Fig. 1141 is oscillatory as well as repulsive. 
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which is again due to the dominant magnetic interaction in the transverse direction. So, with the increase of v the 
magnetic interaction due to the transverse part of the response function plays a crucial role and its form deviates 
from the Lennard- Jones type. This could have important consequences on various binary states in the QGP. 

In QCD the interaction between color charges in various channels is either attractive or repulsive. A quark and an 
antiquark correspond to the sum of irreducible color representations ^: 3(E)3 = 1®8, where the interaction strength 
of the color singlet representation is —16/3 (attractive) whereas that of the color octet channel is 2/3 (repulsive). 
Similarly, a two quark state corresponds to the sum of the irreducible color representations: 3 3 = 3 6, where the 
antisymmetric color triplet is attractive (—8/3) giving rise to possible bound states 45]. The symmetric color sextet 
channel, on the other hand, is repulsive (4/3). Color bound states {i.e. diquarks) of partons at rest have been claimed 
by analyzing lattice data |4a|. The situation is different when partons are in motion. The dipole potentials along the 
direction of propagation and also normal to it have both attractive and repulsive parts, similar to the Lennard- Jones 
form. So, all the attractive channels or the repulsive channels in the static case get inverted due to the two comoving 
partons constituting a dipole in the QGP. This could lead to dissociation of bound states (or to resonance states) as 
well formation of color bound states in the QGP. 

Recent lattice results (sj using maximum entropy method indicate that charmonium states actually persist up to 
2Tc and there are similar evidence for mesonic bound states made of light quarks as well Within our model 
such bound states as well other colored binary states in the QGP can experience different potentials along the dipole 
direction and the direction normal to it. Along the direction of motion binary states which were bound in the QGP 
may become resonance states or dissociate beyond Tc . The temperature up to which they survive need further analysis 
of the bound states properties in detail, which is beyond the scope of this work. Similarly, those colored states which 
were not bound initially in the QGP, may transform into bound states. There are some long distance correlations 
among partons in the QGP, which could indicate the appearance/disappearance of binary states in the QGP beyond 
Tc. On the other hand, in the transverse direction the mesonic states as well as colored bound binary states may 
dissociate for smaller velocity of comoving partons whereas they may remain loosely bound if the comoving partons 
are ultra-relativistic. One immediate consequence of this would be the modification of the transverse momentum, p±, 
dependence of J/ip spectra. Our analysis suggests that a QGP with moving partons beyond Tc does not behave as a 
free gas of partons but a long range correlation may be present. 



III. CONCLUSIONS 



We have investigated the response of the QGP to a fast moving parton within the HTL approximation valid in the 
high temperature limit. For velocities smaller than the phase velocity of the longitudinal plasma mode v — yJiJZc the 
wake potential corresponds to dynamical, anisotropic screening as discussed already before [23|. For larger velocities, 
however, we found, in contrast to previous investigations |23l |24| that Mach cones and Cerenkov radiation appear in 
addition. They have been identified by oscillations (sign flips) in the induced charge density and also in the wake 
potential in the outward flow. This observation could be of phenomenological interest for the present experimental 
programs in RHIC, where various approaches are being followed to flnd Mach cones through the correlation between 
high p± secondary particles in pseudorapidity and azimuthal space. The general structure of the wake potential in the 
direction of motion shows a long range correlation among the particles in colored plasma. The consequences for the 
formation of a possible liquid phase of the QGP |4^J should be investigated. In addition, the wake potential normal 
to the motion of the charged particle causes a transverse flow in the system which might contribute to the creation of 
instabilities |47j and an anomalous viscosity '4^. 

Furthermore, we have investigated the dipole potential by comoving (anti)quarks in the QGP due to the appearance 
of a minimum in the wake potential of a moving color charge, extending the work in Ref . |24| . In Ref. p3 | the dipole 
potential was obtained by averaging the one body potential and the dipole potential was found to be dependent only 
on the electric interaction. In the present work we followed a proper treatment, as depicted in subsec. Ill Fl and the 
dipole potential was found to depend on both electric and magnetic interactions. As discussed both interactions play 
a crucial role in determining the nature of the dipole potential. Based on this we analyzed the various possibilities, 
viz., the appearance/disappearance of bound states in the QGP by discussing the effects of Ampere's law and the 
nature of electric and magnetic forces due to dynamical screening. 



The strength of the color interaction i4Sil can be calculated using SU (3) color group. 
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